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Abstract 

Many-qubit entanglement is crucial for quantum information processing although its exploitation 
is hindered by the detrimental effects of the environment surrounding the many-qubit system. It is 
thus of importance to study the dynamics of general multipartite nonclassical correlation, including 
but not restricted to entanglement, under noise. We did this study for four-qubit GHZ state under 
most common noises in an experiment and found that nonclassical correlation is more robust than 
entanglement except when it is imposed to dephasing channel. Quantum discord presents a sudden 
transition in its dynamics for Pauli-X and Pauli-Y noises as well as Bell-diagonal states interacting 
with dephasing reservoirs and it decays monotonically for Pauli-Z and isotropic noises. 
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I. INTRODUCTION 


fl 


It is aspjmed that eatang.emeat caa be useful iu quautum comput.ug W quautuut etyp- 
togtaphy M aud quantum information processing such as superdense coding H and quantum 

teleportation In practical applications, there are quantum tasks with no entanglement 
but representing quantum advantages over classical counterparts js-?]. This indicates that 
entanglement is not the only feature of quantum correlation and we need to study corre¬ 
lations in a more general concept. This statement is further more correct once it is an 
experimental implementation of quantum computing or quantum information processing. 
In any experiment in the held, entanglement is generally difficult to be achieved, but non- 
classical correlation is comparably less challenging. Due to this fact, and considering the 
additional fact that entanglement in larger Hilbert spaces is not well-known yet, understand¬ 
ing the dynamics of nonclassical correlation, generally, and entanglement, specihcally, is of 
importance specially from practical point of view. 

Experiments are nowadays elaborated such that they can, at least in some physical sys¬ 


tems such as NMR js] and EPR/ENDOR 


0], achieve controls on larger numbers of quits. 


Therefore, it is getting of practical importance to study the characteristics and behaviour of 
larger networks of quits and determine the corresponding dynamics of the correlated states 


under noises that are more common in the physical systems of the stuc 
Quantum discord is one of the measures of nonclassical correlations 


Zl. 


10l |. Its dehnition is 


based on the difference between two classically equivalent forms of mutual information. It has 
been shown that quantum discord is very useful to describe correlations involved in quantum 
systems llj. Nonclassical correlation and entanglement can have in principal different 
behaviours. For example, in some temperature ranges, quantum discord of a two-qubit 
one-dimensional XYZ Heisenberg chain in thermal equilibrium increases with temperature 


while entanglement decreases with temperature and even it goes to zero 


B 


Datta et al. 


showed that entanglement in DQCl is negligibly small, however DQCl includes nonzero 
discord. Thus, the notion of quantum speed up is related to nonclassical correlation rather 
than to entanglement |^. 

Implementations of quantum computing and quantum information processing are ulti¬ 
mate goals for studies and researches in the held thus there have been continuously extensive 
efforts in the context of performing experiments relating to quantum technologies. One im- 
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mediate consideration that is imposed right after starting plans for experimental practices is 
how to deal with decoherence of quantum systems. Decoherence is generally critical for any 
physical system that is determined for experiments, but in the context of quantum proto¬ 
cols this issue becomes one of the most challenging steps. Decoherence of quantum systems 
are due to unavoidable interactions with the environment and it is led to the degradation of 
quantum correlations and limitations on using nonclassical correlation and more importantly 
quantum entanglement in practical applications. 

Studies on quantum channels are usually divided into two categories: Markovian Q,Q 
and non-Markovian [h, 16, 17|. Investigations of nonclassical correlations under Markovian 
environments show that nonclassical correlation is more robust than entanglement. Conse¬ 
quently, implementations of quantum algorithms that rely on nonclassical correlation, and 
not an absolute entanglement, are less fragile 18|. Studies on non-Markovian errors show 


that since there is no occurrence of sudden death of quantum discord in spite of the entangle¬ 
ment sudden death, the description is more practical in the case of nonclassical correlation 

111 - 


Time evolutions of nonclassical correlation and entanglemem 


pled to external environments have been extensively studied 


in bipartite systems cou- 


19|. However, this concept is 


yet not very well resolved for higher dimensions. Quantihcations of nonclassical correlation 
and entanglement for multi-partite systems are not yet generally d known, in addition and 
consequently determinations of their behavior in presence of environmental noises are far 
from being clear. In the case of tripartite systems there have been some results 20|, |23(] but 
for larger spaces the sources are very few 2J-l26| and the existing papers only study the 
entanglement of the states not nonclassical correlations in terms of e.g. quantum discord 
and their behavior under noise. From the previous studies, we have learned that nonclas¬ 
sical correlation and entanglement are preserved in three-qubit GHZ state under noise, but 
correlations are eliminated in the case of a W-state 221. 

In this paper, we study nonclassical correlation and entanglement of four-qubit 
Greenberger-Horne-Zeilinger (GHZ) state that is initially prepared in a pure state described 
as IGHZ 4 ) = ^ (|0000) -|- |1111)). The time evolutions are studied for the IGHZ 4 ) state 
under experimentally common noises, Pauli-X (bit-flip), Pauli-Y (bit-phase-flip), Pauli-Z 
(phase-flip) and isotropic (depolarizing) noises. For each type of noise, the time evolution of 
the state is given by the solution of the master equation in the Lindblad form. The entan- 
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glement evolution of the mixed state is described by using the lower bound for multi-qubit 


concurrence proposed by Li et ah [27|. Nonclassical correlation of the mixed four-qubit 
GHZ state is described by the global quantum discord that is introduced by Rulli et al. [28 1 
which uses a systematic extension of bipartite quantum discord. 

This paper is organized as follows. In the next section, we briefly explain the measures 
that we use for quantifying entanglement and nonclassical correlation of IGHZ 4 ) state. In 
Sec. m we obtain time evolutions of entanglement and nonclassical correlation for the 
initially prepared IGHZ4) state but imposed to Pauli and isotropic noises. In Sec. lYl we 
present the conclusions. 


II. ENTANGLEMENT AND NONCLASSICAL CORRELATION OF IGHZ4) 
STATE 


There is not any general approach in quantifying entanglement for multi-partite states. 
However, as a practical method, one can make use of entanglement measures for bipartite 
states, in which characterizing entanglement is simpler, and extend them to all states by 
convex roof This approach is generally difficult to be performed analytically. Here, in 
order to characterize entanglement of IGHZ4) state, we apply the lower bound for the multi¬ 


qubit concurrence proposed by Li et al. 


27l | which is based on the bipartite concurrences of 


the multi-qubit system corresponding to all possible bipartite cuts of the iV-qubit system. 
For a pure A-qubit state, the concurrence is dehned as follows 


CNim 



( 1 ) 


where pi = Tr(|'0)('0|) is the reduced density matrix of the Tth qubit after tracing out other 
N — 1 qubits. The concurrence for mixed iV-qubit states can be generalized as follows 


Cj^{p) = min E ( 2 ) 

i 

in which the minimum is over the pure-state ensemble that specihes p, namely, p = 
^•Pi|'0j)('0i|. Note, that there is no analytical solution to optimize the concurrence of 
a multi-qubit system except for two-qubit systems and a special case of three-qubit state. 
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Recently, Li et al. 27|] presented a lower bound to this measure which is simpler to calculate 


Cn{p) > T'iv(p) 



N K 




( 3 ) 


This bound contains N bipartite concurrences that correspond to the possible bipartite 
cuts of the multi-qubit system. In order to calculate each bipartite concurrence C'"', one 
needs to evaluate a sum oi K = — 1) terms Ck that are given by 


Cl = max{0, Xl-Xl-Xl- A^}. 


( 4 ) 


Here, A^’s {i = 1..4) are the square roots of the four non-vanishing eigenvalues of the matrix 
PPI in decreasing order and pl = Slp*Sl in which p* is the complex conjugate of p and 
Sp = LI® Lq, k = 1 ,..., K. Moreover, Lq is the generator of the group 50(2) and L^’s are 
the generators of the group S0{2^~^). 

Consider the density matrix of a composed system AB being denoted by p. p^ = Tr b{p) 
and p^ = Tr^(p) are the reduced density matrices for subsystems A and B, respectively. 
Quantum discord for this system, according to the dehnition of Ollivier and Zurek, is given 

by 

D{p) = I{p)-C{p), (5) 

where I{p) = S{p^) + S{p^) — S{p) is the quantum mutual information that is a well-known 
measure of the total correlation and S{p) = —Tr(plog 2 p) is the von-Neumann entropy. In 
addition, 

0(p) = 5(p“)-min{n,}^(p|{n.}) (6) 


is a measure of classical correlations in which the conditional entropy is dehned as 


5(p|{nfc}) = Y,kPkS{pk) where pk = ® ® n^) is the conditional density 

operator, pk = Tr[(J"^ ® nf)p], and the minimum is taken over the set of projective mea¬ 
surements {Hfc}. This form of quantum discord is appropriate for bipartite systems and for 
multi-qubit states we resort to a proper generalization of it. 


There are various approaches to generalize quantum discord for multipartite states 


30 


32|. Here, in order to evaluate quantum discord for a four-qubit system, we apply the 


global quantum discord which is introduced by Rulli et al.. It uses a systematic extension of 


bipartite quantum discord 


28| . The global quantum discord D{paj^...An) for a multipartite 
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state pai...An is then defined as 


N 

D{pA^...A^) = niin{n,}[5(pAi...yi^||<h(pAi...A^)) -'^S{pAj\\<^j{pAj))], 

i=i 

where 5'(pyii,,,Ajv||*^*(PAi...Ajv)) is the relative entropy and 


®(PAi...Ajv) 

hlfc 


i 

y^nfcPAi...Ajvnfc, 




n^/, 

An' 


(7) 


( 8 ) 

(9) 

( 10 ) 


and k denotes the index string (ji...jAr). 

In order to define local projective measnrements, we select a set of von-Neumann mea¬ 
surements as 

^1 ^ / cos^d) cos(|) sin(|) 

\^e-*'^7cos(|)sin(^) sin^d) 

and 

^2 ( sin^(y) cos(^) sin(^) 

^ cos(^) sin(^) cos^(^) 

where 6j G [0,7r) and (j)j G [0,27r). By varying the angles 6j and for j = 1..4, one can 
find the measurement basis that minimizes Eq. ra¬ 




in. TIME EVOLUTION OE IGHZ4) STATE UNDER NOISE 


Time evolution of a quantum system under noise is given by the master equation in the 


Lindblad form 


33| 


dp 

dt 


^[Hs,p\ + [Ci^apC 


'la 







(13) 


where Ci^a = s/i^^crl is the Lindblad operator that describes the noise and acts on the 
ith qubit. is the decoherence rate and al are the Pauli spin matrices for the ith qubit 
with a = X for Pauli-X, a = y for Pauli-Y and a = z for Pauli-Z noise. Also, is the 
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Hamiltonian of the system. Solutions to the above equation under various noises for four- 


qubit GHZ state is recently presented in Ref. 


34j |. The idea is to write the density matrix 


for inhnitesimal time interval t = 6t hj using the Lindblad equation as 

1 


p{6t) = p(0) 


5 ; (£,,„p(0)4„) - - 


5t, 


(14) 


where p(0) = |GHZ 4 )(GHZ 4 |. Then, by using a proper ansatz, the solutions can be obtained 


for all time 


34|. 


Under different noises, if each qubit interacts locally with the environment then the 
dynamics of nonclassical correlation and entanglement can be extracted from the above 
formulation, specihcally for an initially prepared four quit GHZ state. 

A. Pauli-X noise 


The density matrix of four quit GHZ state under the Pauli-X error is the solution of the 


Lindblad equation [3^, as follows 




a 0 
0 /3 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 /3 

a 0 


0 0 0 0 0 

0 0 0 0 0 

/3 0 0 0 0 

0 7 0 0 0 

0 0 /3 0 0 

0 0 0 7 0 

0 0 0 0 7 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 7 

0 0 0 7 0 

0 0 /3 0 0 

0 7 0 0 0 

/3 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 7 
/3 /3 0 
0 0 0 
0 0 7 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 


0 0 0 0 0 a 

0 0 0 0 /3 0 

0 0 0 /3 0 

0 0 7 0 0 

0 /3 0 0 0 

7 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

7 0 0 0 0 

0 /3 0 0 0 

0 0 7 0 0 

0 0 0 /3 0 0 

0 0 0 0 /3 0 

0 0 0 0 0 a 


\ 


(16) 


where 


a ■■ 

f^- 

7 = 


16 

J_ 

16 

1^(1 


(1 + + e-®''*), 

(1 - e-®''*), 

■ + e"®^*), 


(16) 


and K is the decoherence parameter. For this mixed state, the lower bound Eq. ([3]) to the 


four-qubit concurrence gives 


r(p^) = max <j 0, ^ (e"®"* 


-|- 6e 


— 4Act 


3 ) 


( 17 ) 


which is an exponentially decaying function of time. The plot of this lower bound (Fig. [T]) 
shows that the entanglement of the four-qubit GHZ state vanishes after some hnite time for 
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FIG. 1: The lower bound for the four-qubit concurrence for an initial IGHZ 4 ) state, transmitted 
through Pauli-X, Pauli-Y, Pauli-Z and isotropic channels as fnnction of Kt. 


the bit-flip channel. However, using the positive partial transpose criteria 29|, we hnd it 
out that the density matrix of ffT^ is separable only for t ^ 00. This indicates that the 
introduced lower bound does not display proper long-time entanglement for this state. 

In order to obtain nonclassical correlation, we need to evaluate Eq. ([7]) by hrst Ending 
the reduced density matrix of ith subsystem which is obtained by tracing out other three 
qubits that reads 

p\ = p\ ^ P\ ^ dA4 = (a + 4/3 -b 37 ) 1. (18) 


Since p\. is proportional to the unity operator, is unafiected under 

any projective measurements, therefore, = Pa^ and the second term in Eq. ([7]) 

vanishes. To evaluate the first term, we need to find the measurement bases that minimize 
quantum discord. Calculations show that quantum discord suddenly changes during the 
dynamics of the system, namely, the measurement bases change from {9 = 0) to 
{6 = Ti 12), so we find 

^($(p^)) = 4 - ^^^ log,(l - 

-^[(1 + + e-®"') log2(l + -f- e"®"*) 

8 

+3(1 - log2(l - 26"^^* + (19) 


for 0 < Kt < 0.137 and 

^(«l>(p^)) = 3, 


( 20 ) 
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for Kt > 0.137. Also, 


n _ - 8 Kt\ 

S{p^) = 3 - ^^ log 2 (l - 

-^[(1 + 6e-^"* + e-®"*) log2(l + 6e-^"* + e"®"*) 

8 

+3(1 - + e-®"*) log 2 (l - 26"^"* + e"®"*)]. 


Finally, we obtain 


D{p^) = 1, 


for 0 < Kf < 0.137 and 

Dip^) = - -^-^log2(l-e-«-*) 

+ ^[(1 + + e-®"') log2(l + 6e-^"* + e"®"*) 

8 

+ 3(1 - 26-“ + 6-®“') log 2 (l - 26-““' + 6-®“')], 


( 21 ) 


( 22 ) 


(23) 


for nt > 0.137. Therefore, the quantum discord is constant for 0 < Kt < 0.137, however, 
it experiences a sudden transition in its dynamics at nt = 0.137 and monotonically decays 


for Kt > 0.137. This peculiar behavior of quantum discord is known as t 
of the quantum discord which is also observed for two qubit systems 
universal freezing of quantum correlations 37|. 


as t 

le 

35, 

CO 


le sudden change 


36l |. linked to the 


B. Pauli-Y noise 


Now, we investigate the time evolution of nonclassical correlation and entanglement for 
IGHZ 4 ) state under Pauli-Y error. The density matrix for this system is given by 341 




P = 


/ 

a 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

/ 

0 

P 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-p 

0 


0 

0 

p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-p 

0 

0 


0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

-7 

0 

0 

0 


0 

0 

0 

0 

p 

0 

0 

0 

0 

0 

0 

-p 

0 

0 

0 

0 


0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

7 

0 

0 

7 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

p 

-p 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

-p 

p 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

7 

0 

0 

7 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 


0 

0 

0 

0 

-p 

0 

0 

0 

0 

0 

0 

p 

0 

0 

0 

0 


0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 


0 

0 

-p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

p 

0 

0 

1 

0 

-p 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

p 

0 

\ 

a 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 


(24) 


/ 

For this case, the lower bound to the concurrence is similar to Eq. (UBi i-e-, = 

max |o, ^ ~ 3)|. Also, similar to Pauli-X noise, we observe a sudden change 
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in quantum discord at nt = 0.137. Indeed, for nt < 0.137 the optimal local measurements 
are given by 11^^. = |i)(i| and since ^j{pA.) = Pa. and — S{p^) = 1, we find 

D(y) = 1 . (25) 

But, for Kt > 0.137, ay is the basis which minimizes quantum discord, so we have S'(<h(p'^)) = 
3 and S{p^) = S{p^). Therefore, D{p^) coincides with Eq. fl23D . 


C. Pauli-Z noise 


For the next case, consider Pauli-Z error where its corresponding density matrix takes 
the following form 341 




+ |1)®^(1|®^) + -e-®"* (|0)®^(1|®^ + |1)®^(0|®^) , 


(26) 


that results in the lower bound 

t{p^) = ( 27 ) 

As it is shown in Fig. [H nonclassical correlation and entanglement of this state decrease 
more slowly with respect to the previous noise models. For this case, the reduced density 
matrices become 

Pa, = Pm = Pm = Pm = \ ( 28 ) 

which results in <1)^ {p \.) = pA ^, therefore the second term of Eq. ([7]) vanishes. The optimiza¬ 
tion procedure of 9 and 0 for each qubit shows that the optimized projective measurements 
are the bases of a^- Then we hnd S{^{p^)) = 1 and the von-Neumann entropy becomes 

s(p") = 1 - ^(1 - log 2 (l - e-®”') + 01 + e-*'*) logj(l + e-®”'). (29) 

Therefore, we have 

D(P^) = 01 - e-*“‘) logj(l - e-M + (1 + e-*“‘) logj(l + e-*”'). (30) 

This result shows that the Pauli-Z error causes dissipative behavior on nonclassical correla¬ 
tion for IGHZ 4 ) state. Decoherence behavior in terms of nt is depicted in Fig. | 2 l 
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D. Isotropic noise 


To this end, consider IGHZ 4 ) state under isotropic noise. It should be noted that for this 
case, the Lindblad equation depends on all three Pauli operators (see Ref. 29). Thus, the 
density matrix as the solution of the master equation is described by 34l | 


/ 

a+ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a_ 

\ 

0 

/3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

/3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

/3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

/3 

0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

/3 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

/3 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

p 

0 

0 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

p 

0 


\ 

a— 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

«+ 

/ 


where 


' a+ = ^ (^1 + be-®''* + , 

7 = ^ (^1 - , 


(32) 


Now, the lower bound ([3]) to the concurrence reads 



max 


0. # (9. 


— lQK,t 


+ 6e 


— Stvt 



(33) 


This lower bound shows strong dissipation of entanglement in presence of isotropic error. 
It should be noted that for this density matrix, similar to the Pauli-X and Pauli-Y errors, 
according to the positive partial trace criteria, expression fl3^ is only valid for t < to. 

The reduced density matrices become 


Pai =p52 =Pa3 =^^4 = (« + 4/^ + 37)1- (34) 

So, Eq. ([7]) reduces to 

D{Pm...aJ = min{n,}>^(p!4i...A^II*^’(PAi...A^))- (35) 
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Calculations show that the optimized local projective measurements are = |i)(i|. Thus, 
we have 

(-i _ -lQKt\ 

5(p°) = 4-^-^ 

--(1 - log 2 (l - 

8 

- — (1 + 6 e-®"* - logofl + de"®"* - 

16 

+ —(1 + 6 e-®^' + log 2 (l + (36) 

16 

and 

^(«l>(pD)) = 4-^-^-^log 2 (l-e -^®-0 

--(1 - 2 e-®'^‘ + log 2 (l - 2 e-®''' + 

8 

-^(1 + 6 e-®"* + log 2 (l + 6 e-®"* + (37) 

8 

By inserting Eqs. fl36|) and flTTI) in Eq. fl35|l then 

D(pD) = -l(l + 6 e-®''* + e-^®''*)log 2 (l + 6 e-®^* + e-i®''*) 

8 

+—(1 + 6e-®''* - 7e-^®'^*) log2(l + de"®'^' - 76-^®''*) 

16 

+—(1 + de-®''* + 9e-^®'^*) log 2 (l + de-®'^* + 9e-^®'^')]. (38) 

16 

The time evolution of quantum discord for this case is plotted in Fig. O As the hgure shows, 
the isotropic error results in a more decrease in quantum discord in comparison to Pauli-Z 
channel. 

IV. CONCLUSIONS 

In this paper, we investigated the dynamics of nonclassical correlation and entanglement 
for the initial IGHZ 4 ) state under interactions with several independent Markovian environ¬ 
ments. These noises include Pauli-X, Pauli-Y, Pauli-Z and the isotropic errors. In order 
to describe the time evolution of entanglement, we utilized the lower bound ([ 3 ]) as an ap¬ 
proximation for the convex roof of the four-qubit concurrence. The results showed that 
entanglement is more robust in transmission under the Pauli-Z error compared to Pauli-X, 
Pauli-Y and isotropic errors. For Pauli-X, Pauli-Y, and isotropic errors the lower bound 
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0.0 0.2 0.4 0.6 0.8 1.0 


FIG. 2: Quantum discord for the four-qubit system with the initial IGHZ 4 ) state transmitted under 
Pauli-X, Pauli-Y, Pauli-Z, and isotropic errors as a function of nt. 

vanishes in a finite time (sudden death phenomenon). However, from the positive partial 


transpose separability criteria 


29| . we conclude that these states are separable only for 


f —)■ cxD. Thus, the lower bound only describes the behavior of the entanglement for t < to 
where the lower bound vanishes at to- This is due to the fact that the accuracy of the lower 
bound over convex roof seems to depend on the rank of the corresponding density matrix 


2l| . Indeed, the numerical calculation of the convex roof is a quite difficult task which re¬ 


quires an optimization over too many (proportional to where r is the rank of the density 
matrix) free parameters 211 . 

In order to obtain the time evolution of nonclassical correlation, we used the global 
quantum discord as a measure to calculate nonclassical correlations in multipartite states. 
We found that quantum discord is not affected by Pauli-X and Pauli-Y noises for nt < 0.137 
whereas it decreases monotonically for nt > 0.137. This transition of quantum discord is due 
to the transition of measurement bases which minimize the quantum discord. Also, quantum 
discord decays in the presence of Pauli-Z and isotropic noises in terms of the dimensionless 
scaled time nt. In comparison, quantum discord oflGHZa) state is more robust than IGHZ 4 ) 
against decoherence for all noises under study [38|. Moreover, except the case of Pauli-Z 
error, quantum discord is more robust than entanglement for the initially prepared IGHZ 4 ) 
state. 
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